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I- SUMMARY OF ACCOMPLISHMENTS 


Assembly of data and hydrologic and physiographic 
descriptions of the Soviet and U.S. test sites- 

The following documents nave been prepared for presentation 
to the NASA technical monitor: 

a. A geographic description of the Fergana Valley 
test site was sent to the NASA Technical Monitor in 
April 1977- This description will be updated, including 
information from the Soviet description, for the final 
report. 

b- A geographical description of the southern 
Sierra Nevada/San Joaguin Valley test site was sent to 
the NASA Technical Monitor in June 1977- 

c- Numerous field measurements of solar 
radiation and snow reflectance taken during the 
1976-77 snow season have been tabulated and stored as 
a computer text file- Printed copies of this data are 
available upon reguest and will be included with the 
final report- 

d» Field measure in ants of air temperature, vapor 
pressure, surface temperature, and incoming 
atmospheric radiation, taken during the 1976-77 snow 
season, have also been tabulated and stored as a 
computer text file- Printed copies of this data are 
available upon reguest and will be included in the 
final report- 


2. Development of energy-balance snowmelt model for 
large mountainous area- 

a. Development of wavelength-dependent topographic 
solar radiation algorithms- 

i- Method of calculating absorption and 
scattering coefficients and portion of radiation 
scattered forward has been completed- 

ii. These coefficients and terrain 
information are used to produce maps of solar 
radiation. At present the algorithm does not 
include non-Lambertian reflectance from adjacent 
slopes, but this modification is under development. 

iii. Numerous (over 400) field measurements 
of incoming and reflected radiation have bean made, 
for calibration and verification of algorithms- 
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b, Determimation of snow reflectance from Landsat- 

i. We have adapted DIRS software to our 
computer system and have identified suitable ground 
control points within study area- 

ii. h program to map space radiance values 
corresponding to incoming solar radiarion has been 
developed. 

iii. We have necessary field data to apply 
standard atmospheric corrections- 

c, Ilemperature and humidity estiraates. 

i. An algorithm for interpolating temperature 
and vapor pressure over a topographic surface from 
a few discontinuous measurements has been 
developed, 

ii. fieasurements of temperature and vapor 
pressure have bean made at various locations, and 
elevations in the study area, at times coincident 
with Landsat passes throughout the snow season. 

Data from the Bishop HSO has been combined with 
these measurements for utilization of the above 
mentioned algorithm. 

d- Development of topographic longwave radiation 
algorithms. 

i, a method for calculating incoming longwave 
radiation under clear skies from near surface 
temperature and vapor pressure has been developed. 

ii. Field measurements of incoming longwave 
radiation have been taken using an Sppley 
pyrgeometer for verification and calibration of 
algorithms- 

e. Wind interpolations - We have not modeled 'wind 
over 3- diraensicnal terrain. But we have begun to devise 
some novel methods of indirectly measuring wind, using 
thermal imagery. 


3. Coupling of smowmelt model to a flow model- 

a. We have examined the literature on frost 
penetration ‘and its effect on infiltration rates. 
During October we plan to install thermistor atfays at 
several locations in the study area, so that we can 
monitor soil and snow temperatures during the winter. 

b. We have begun to develop, a geocoded hydrologic 
information system for the study araa. 


c- We have examined literature on flow through 
snow, and have begun to formulate methods of routing 
meltwater to the streams,. 


4, Review of water balance and hydrologic literature 
to determine potential contributions of remote sensing. 

a- We have assembled a set of about 500 pertinent, 
recent references, and have developed a computerized 
literature retrieval system around rhera. 

b. Future references can be easily included in our 
system, and it will shortly be expanded greatly with our 
collection of soil moisture literature- 

c. During the remainder of the first year of the 
grant period we will expand our review to include other 
potential applications of remote sensing to hydrologic 
modeling. 


II. DISCUSSION OF TASKS UNDEEWAY 


A. AN ENERGY BALANCE SNOW'-ISLT EODEL FOR RUGGED TEPSAIN 

One of the major objectives of the project is the 
extension of an energy balance snowmelt model over rugged 
terrain. In the initial season we have primarily 
concentrated on radiation inputs into the snowmelt model, 
and we have made progress in deriving methods of calculating 
solar radiation and longwave radiation, and in determining 
wavelength-dependent reflectance from Landsat data. We have 
also made a very large number of field measurements of solar 
and longwave radiation, temperature, humidity, wind speed, 
and snow temperature and d.snsities; those measurements will 
be used to calibrate and/or verify our computations- By the 
end of the first year of the grant hope to have the 
radiation algorithms solved- These studios are discussed in 
detail below. 

We have also worked on ways to derive wind velocities 
over a large area, and ways to couple a snowmelt model to a 
flow model. In future field seasons «e hope no work on 
these probleras- 


1- Development and Testing of v/avelength-Dependent 
Topographic Solar Radiation Algorithms 

Many slope radiation models have been developed, but 
typically they either do not include horizon information, do 
not calculate reflection from adjacent slopes, or do not 
subdivids incoming radiation into wavelength bands. In a 
mountainous, snow-coverai area, all of these refinements are 
essential. 

Here we have developed a program which use;! the NCIC 
Digital Terrain Tapes to calculate slope, exposure, and 
horizon information for every point on a grid (Dozier and 
Outcalt, in preparation). 

To calculate the solar radiation distribution on such a 
surface, it is necessary to first calculate the beam and 
diffuse radiation at the range of altitudes in the region. 
Generally the approach to this problem has been to model the 
incoming radiation using estimates of the atmospheric 
aerosol distribution, but these methods are computationally 
expensive, and aerosol distributions are typically not 
known. 

Our approach is to calculate the wavelength- dependent 
absorption (hx) and scattering coefficients and the 

portion of the radiation scattered forward (S^) from a set 
of three global radiation measurements with an Exotech 
Landsat Ground-Truth Radiometer- The method used is to take 
three measurements, which differ either in orientation or in 
optical air mass (which varies with solar zenith angle Z) . 


The team radiation on a 


surface normal to sun is: 


^ ^0 I - ■*‘^1 “'s 5 J 

where S-^q = solar constant in band X 

mg = relative absorbing air mass, either 
ozone or water vapor, depending on X 
mg = relative scattering air mass. 

The beam radiation on a horizontal surface is; 


” ^Xi. ^ 

Global radiation on an unobstructed horizontal surface is 
(Kondratyev, 1969, p- 463) : 


Gxh = i2„-_I!!s.ULSxx).-.^KElzi£xIfe.l-.£SS_2J, 

2 (1 - " 2^*XD5 ®^PC^A^X ^® s “ ^ 

where A;;y^p is the diffuse reflectance in wavelength band X, 
The diffuse radiation on a horizontal surface is 

thus; 


“ ^Xh ”^Xh 

The global radiation at a point on a slope, which is the 
value measured in the field, is the sum of the beam, 
diffuse, and reflected components: 


'^Xsl ■= SXsT + D;)^s-| + B.xs-\ 

where - 



Sxj_ cos Z V 



%sl " ^Xh " 
^X1 


sin 9” 


(1 



where Z* = 
Vf = 

^XDB = 

?-XB = 

H = 
sin 0" = 
S' = 
Z'» = 
<p = 
<p. ^ 


solar zenith angle with respect to slope 
view factor, portion of sky seen, 

= cosz (0-5*[S+ (90 - H) J) 

portion of beam radiation which is diffusely 
reflected 

portion of beam radiation which is specularly 
reflected 

mean horizon angle from zenith 

Cos S' cos Z" - cos Z" sin Z" cos (fp-ip') 

slope with respect to opposing slope 

solar zenith angle on opposing slope 

azimuth of sun 

azimuth of slope 


The formulation for sin 9" is from paitridge and Platt 
(1976, p,128) - 


..Nl , L. \ .1 1 j I I: ^ f 
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Three measurements of the global radiation, taken either 
at one time with different orientations, or taken at 
different times and hence different 2, lead to three non- 
linear eguations with three unknowns, These are solved using 
Brown’s (1969) method, vJe will use our extensive radiation 
measurements to evaluate the effectiveness of this algorithm. 


2. Shortwave Albedo of an Extensive Snow Cover Derived 
from Landsat Data 


Data inputs to the model consist of the following; 

a. Field shortwave radiation measurements: We use 

Exotech model 100-A "Landsat Ground-Truth" radiometers, with 
swltchable spectral bandpasses identical to MSS bands 4 
through 7- The current data base contains measurements taken 
at intervals of approximately two weeks from 12 February to 
31 Kay 1977, coinciding with Landsat overpasses of the study 
area. The measurement sequence is described in the previous 
section. In addition we make a reading of the sky radiation 
with a 15° field-of- view lens, 

b, Landsat computer compatible tapes (CCT's).; scenes 
coinciding with field measurements are a minimum reguirement. 
Ideal coverage requires 3 scenes for each overpass interval 
of 2 days. Images acquired in the high-gain mode are 
essentially useless owing to detector saturation over snow, 

c- Digital Terrain Tapes; topography of the study area 
(digitized at 63,5-raeter intervals) is obtained on tape from 
the national Cartographic Information Center, 

d. Topographic maps: at 1:62, 500 (15') scale, they 

provide complete coverage of the study area, from which 
geodetic coordinates of any locatable feature or measurement 
site may be manually extracted. 

e. Ephemeral data; a computer program simulates a 
specific ephemeris for use in sun angle calculations. 

Certain primary transformations must be performed upon 
each data input. These are; 

a. Field shortwave radiation ineasurements are converted 
from analog voltage readings to irradiances (W m"^) aud 
spectral irradiances (W m-^^ nin“*-) , by program RADFRNT, This 
same program organizes the measurements into an indexed file 
from which they may be accessed according to several sorting 
parameters. Albedos (beam and global) are calculated by 
ratioing incoming and outgoing values in subsequent programs. 

b. Landsat digj.tal voltage readings (from CCT's) ate 
converted to space upward radiance values {H(space)) 

(W m ~2 sr"'!-) by the following algorithm. 

R{space) = [ Vc-v: (R (max) -R (min) ) / (Ks*Ko) j + R (min) 


1 


NL 1 . r 1 1 j i 
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Vc is the calibrated digital voltage for each pixel for 
each of 4 bands- R [max) and R (min) are band- dependent 
saturation and threshold radiance values. Ks is a 
sun-correction constant. Ko is 127 for compressed data from 
bands 4 through 6, 63 otherwise. In practice this algorithm 
is implemented as a table lookup by program RIDVAL, 

Saturated values are set to 0 , and all xero values are 
subsequently ignored- RADVAI also reverses the row-order of 
the image so that the origin is in the southwest corner, 
conforming to the project standard. Pixel coordinates ^ 

(sample and line) of selected features are manually extracted 
from transparencies of the Landsat scene, and verified from 
computer shade prints produced by program COMPOS. "These 
surface features are evaluated for possible use as ground 
control points {GCP*s) in subsequent image rectification. 

c- Digital terrain Tapes are processed by program 
DTTOEG, which i) resamples the topographic grid to a standard 
interval (a 1 00 m interval has been selected as an optimal 
compromise between resolution and date storage requirements) ; 
ii) transposes the grid to a west-to-east major order, 
beginning at the south-west corner. This resampled and 
reordered terrain data serves as the base grid to whicn all 
other spatial data is registered. In addition, mean slope 
and aspect are computed and stored for each grid cell. 

d. Topographic maps are used to precisely locate field ; 

measurement sites and those image grid control points for 
which reliable image coordinates have been found. 

Once primary data transformations are completed, the 
following seccndacy transformations are required to complete 
the database. 

a. Relevant portions of the Landsat imagary axe ; 

rectified to a UTM projection and are resampled to the 

standard grid interval of 100 m- This is accoinpiished by 
various programs in the DIRS package, originally written at 
NASA/Goddard SFC, The rectification uses original ground 
control data or is accomplished by registering new imagery to 
a previously rectified image of the same area. 

b. The Digital Terrain data is combined with ephemeral 
data to generate solar -zenith angles for each point on the 
grid. These are computed relative to absolute and slope 
zeniths, for the date and time of each image- The 
slope-zenith angle is combined with ephemeris data and the 
solar constant to yield a spectral value for incoming 

shortwave radiation ignoring the atmosphere, ; 

c- Field shortwave mGasareirients and the corresponding i 

computed maxima are input to an atmospheric model described 
elsewhere in this paper, from which empirical spectral t 

coefficients of absorption, scattering, and fractional 
forward scattering are in turn used to adjust the computed ;i 

incoming shortwave radiation throughout the grid, i 
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d. The same atmospheric attenuation factors are now 
applied to the rectified Landsat radiance values, this time 
as compensation for the attenuation of upwelling radiance. 

We now have a grid consisting of corrected upwelling radiance 
and corrected surface irradiance, computed for each grid 
cell, 

e. The albedo calculated above assumes the surface to 
be an isotropic (Lambertian) reflector. In fact, snow is an 
anisotropic reflector with a specular component of reflection 
that becomes quite pronounced at small angles of incidence 
(large slope~ normal zenith angles). The slope and aspect 
data contained in the grid, combined with ephemeral 
calculations, make it possilbe to introduce a correction to 
the computed albedo of each cell, based on sun angle, mean 
slope and aspect of the cell, and look-angle of the KSS, The 
precise nature of this function has not yet been determined; 
its dependehce upon the age of the snowpack may require an 
empirically-determined relationship. Also required is a 
precise determination of which grid cells actually represent 
unobstructed snow cover, for which an angularly dependent 
albedo correction would indeed be approporiate. 

Of the steps outlined above, the following have been 
accomplished: 

i- collection of sufficient raw data to begin 
a meaningful analysis, 

ii, conversion and indexing of field 
shortwave measurements, 

iii, resampling and regridding of digital 
terrain data, 

iv, extraction of geodetic coordinates of 
image control points and field measurement sites, 

V. calculation of ephemeral data (solar 
zenith angle, solar constant) for entire grid at 
■ measurement and/or imaging intervals, 

vi, rectification and resampling of a Landsat 
subiraage, 

vii. calculation of slope, aspect, and 
horizon profiles for each grid cell. 

Of the remaining steps in the model, the following are 
currently under development (i-e,, projected completion 
within one month): 

i, calculation of "clear-sky" shortwave 
irradiance ever the grid, 

li, refinement of the atmospheric model and 
generation of attenuation coefficients. 


are : 


iii- registration of the rectified Lan3sat 
data to the base (terrain) grid- 

The Biajor steps in the model remaining xo be i>>ipleinented 


i. calibration of the Landsat and surface 
radiance values with the atmospheric factcjs, and 
subsequent albedo calculation, 

ii- development of a workable angular 
correction for snow albedo- 


3- Development and Testing of Longwave Radiation 

Algorithms 

The model under development is structured as follows^ 

a- A digitized grid of elevations is formed for the 
study area from. Digital Terrain Tapes (DTT* 3} using the 
program DTTORG- The grid cosolution is variable and the 
optimum resolution is determined for each model application. 
Mean horizon angles and thermal view factors are calculated 
from terrain data for each grid point by rouxineHORTAB 
(Dozier and Outcalt, in preparation) : 

Vf = cos2[ (90-H) /2] 

where; Vf = thermal view factor 

H = mean horizon angle from grid point 

b,' Field data and data from Bishop WSO are combined by 
routine LNGDAT- Measurements of air temperature, wet bulb 
temperature, air pressure, and elevation are combined 
sequentially with field measurements from a particular model 
run period. These are used as a reference for calculating 
air pressure over the study area and for the establishment of 
temperature and vapor pressure lapse rates- Pressure is 
calculated for each grid point using Bishop station pressure 
and a standard temperature lapse rate (>-7°K km“^ is typical 
of the study area). Note that this equation is reasonably 
insensitive to minor variations in the temperature lapse rate 
due to the extreme variation of pressure with elevation 
(Kantor and Cole, 1 965);: 

Pa = Po ^ exp[ (-g=^=m) / (3*r ) * In ((To + r* (z-zo) ) /To} ] 

where: Pa = pressure at elevation z 

Po = pressure at Bishop 
g = acceleration of gravity 
m = molecular weight of dry air 
S - molecular constant for an ideal gas 
r = temperature lapse rate 
z = elevation of grid point 
zo = elevation at Bishop WSO 
To = temperature at Bishop WSO 
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c. Temperature^ and vapor pressure lapse rates are 
established for the model run period over the study area from 
field measurement data- Wet bulb depression (Ta-Tw) is 
converted to vapor pressure from station pressure and the 
psychrometric equation (3yers, 1974); 

e (a) = e (s) - [ (Cp*Pa) / (, 622*L) ] * (Ta-Tw) 

where: e(a) = vapor pressure 

e(s) = saturation vapor pressure 

Cp = specific heat of dry air at constant 
pressure 

1 = latent heat of vaporization 
Ta = air temperature at z 

Temperature is converted to potential temperature to simplify 
computations (C ramer , 1 972) : 

9 = Ta * (1000/Ba)«>!<[R/(m*Cp) J 

Field measurements of temperature and vapor pressure, which 
are irregular both spatially and temporally, are expanded to 
a regular grid by program N13HGHD using the following 
procedures: 

i. Gravity model forms an irregular grid from 
X/Y relationship of elevation and time- 

ii. Laplacian transform smooths the irregular 

grid (Note: procedures a) and b) were developed by 

S. Outcalt) 

iii- A bicubic spline interpolation resamples 
the irregular grid at specified time and elevation 
intervals forming a regular grid (de Boor, 1962; 

Rogers et al, 1976)- 

These grids of [ (ae/St) by (aa/az) ] and [ (ae{a)/at) by 
(ae {a)/az) J are used to determine 9 and e (a) for each point 
in the elevation grid- Potential temperature is converted 
back to air temperature for each grid point by: 

Ta = 6/ (1000/Pa) **[?./ (m*Cp) ] 


d„ Program INFLUX calculates the following: 

i- Effective clear sky emmissivity for every 
grid point using Brutsaert's (1975) techuigue 
corrected for diminished pressure and optical path 
length: 

€ (a) = 1-24[e(a)/Ta]V‘^ 


ii- Incoming longwave radiation for every grid 
point using output from the above eguations: 

I(i) = [€ (a)^ * [1-VfJ + J * vf 

where: a = Stefan-3oltzraann constant 

Ts = surface temperature 

Surface temperature is approximated from an interpolated 
mean of field measurements- 


e- The grid of incoming longwave radiation flux is: 

i. Written to disk for input to energy balance 
snowmelt model, 

ii. Plotted as a contour map over the study 

area, 

iii. Plotted as a three dimensional surface 
over the study area- 

The model will be calibrated and tested for accuracy by 
the following procedures: 

a. Bishop WSO data will be used to calculate pressure 
over the study area and to initialize the model for Ta, e (a ) , 
and Pa, 

b. Field measurements of la and e (a) will be used to 
determine accurate temperature and vapor pressure lapse rates 
for determining effective atmospheric eramissivity, 

c- Pyrgeometer measurements of incoming longwave 
radiation will be used: 

i- To run a sensitivity analysis of input 
parameters Ta, e (a) , their lapse rates, and Pa, the 
grid resclution, and the time interval, to improve 
the model efficiency and the cost-effectiveness of 
field data collection, 

ii. To make a first order determination of the 
effect of cloud and forest radiation on incoming 
longwave radiation- 

iii. To evaluate model accuracy statistically. 


4 - Possible hoproaches to Wind Interpolation 

One of the knu'^ty problems which to coiiuinuas to- face us 
is the interpolation of wind over complicatad three- 
dimensional terrain. While flow over simplified, two- 
dimensional topography has been simulated, as has flow at 
higher elevations in the atmosphere, calculation of surface 
turbulent exchange remains a serious problem. During the 


next seasons, we intend to try to utilize thermal data from 
Landsat-C, HCHH, and Tiros JI to investigate the wind field* 
Our approach will be to invert the snow surface temperature 
model- Specifically, during non-melting conditions, if wind 
speed, roughness length, radiation etcS are known, a snow 
temperature profile is one of the outputs of Anderson's 
(1976) model- The existence of thermal satellite data would 
allow us to measure surface temperature, and solve for the 
turbulent transfer coefficient, which is related to wind 
speed- 

In order to properly carry out this experiment, we would 
need thermal data 24 hours apart for the same area- For 
Landsat-C or HCHM data we could use the sidelap zone, or we 
could use Tiros N data directly. Landsat-C would provade the 
best resolution, but the temperature response of the 
instrument is not as precise- HCMM data may not be readily 
available. In any case we would carry out simultaneous 
ground measurements, to filter out atmospheric effects- 

Once several spatial measurements under different 
meteorological conditions were obtained, we could estimate 
area- wide wind distributions for these conditions- 
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B- ecu FLING SNOWMELT KODEL TO A PLOW MODEL 


1^ Review of Literature on Frcet penetration and its 
Hydrologic Effects 

We have reviewed the general properties and processes of 
seasonally frozen ground emphasizny relatively coarse soils, 
;-lnce preliminary research of the mountainous portion of the 
test area indicates the predominance of rocky sandy soils 
(Anderson, 19 46, Storie, 1 957)^ Heat and moisture exchange 
between the snowpack and ground can influence the timing of 
snowmelt. Conditions on the surface and in the ground can 
also influence the timing and rates of runoff. 

The presence of bonded concrete frost represents and 
impermeable layer at the depth of ice formation (Dingman, 
1975), Although concrete frost is probably rare in the 
Sierra, it may occur in years when initial snowfall is late. 
Other types of frost may increase infiltration capacity or 
leave it unaltered (Trimble, et al- , 1958, Schuram and Lusby, 
1963). Outcalt (1971, 1976) has developed computer 
simulation of the formation of frost with varying morphology, 
and Harlan (1971, 1973) presents a coupled heat-fluid 
transport model which descrihres the process of infiltration 
into a freezing or frozen soil with given surface boundary 
conditions. The process of infiltration is also nuraericallly 
described in static phases by a numerical model developed by 
Alexeev et al. (1972) who also present a simplified dynamic 
model shown to give good results in sandy soils- The models 
by Harlan (1971 , 1973) and Alexeev et al, ( 1 972) , relatively 
recent and theoretically based, lack extensive empirical 
validation. These models, particularly that of Harlan (1973) 
may be suitable for insertion as subroutines in an overall 
surface energy budget simulator {e,g, Outcalt and Carlson, 
1975) . 

In order to test the ability of these various models to 
calculate soil temperatures and frost penetration, we plan to 
install arrays of thermistors and soil moisture blocks at a 
few locations within the study area this fall. These will 
enable us to rapidly measure temperature and moisture 
gradients within the soil (and snow) when we visit the sites. 


2. Development of Geocoded Hydrologic Information 
System 

We have begun to assemble a geocoded hydrologic 
information system for the study area- The drainage basins 
have been cla5!isified into a hierarchy based upon drainage 
area, and we are' digitizing the boundaries. Once the 
boundaries are cligitized, we will be able to Create a file 
for a given basin which contains a variety of data for every 
point in the basin. Static data would include topographic 
information such as slope, exposure, horizons, elevation, and 
vegetative cover. Temporal data would include solar 


radiation for a giver data, Landsat radianct>, ^aflectaticij 
and whether snow-covered or not. 

We have not yet agreed on the proper format for this 
information system, and plan to delay doing so until after a 
NASA-funded conference on geobase information system formats^ 
to be held in Santa Barbara this month, is concluded. 


3, Development of Flow Model for Snowmelt 

In the Southern Sierra Nevada almost all of the runoff 
is from snowmelt. Under extreme conditions, such as occurred 
during this last winter, snowmelt models based on regression 
analysis fail to adequately predict snowmelt runoff, Energy 
balance snowmelt models can predict runoff accurately under 
extreme (drought or flood) conditions, but to be used 
operationally to improve development or regulation of water 
resources, they roust be coupled to existing or modified 
streamflow models. 

Thus the energy balance snowraelt model must be coupled 
to a strearaflcw model. Once snow has begun to melt the 
following must be evaluated: 1) infiltration and runoff 

beneath the snowpack; 2) moisture flow in saturated soil 
below the snow line; 3) soil vapor losses- In addition 
problems of differences in the resolution between the 
snowmelt model and the flow model must be overcome. We see 
our primary problem in the estimation of flow from the snow 
to the nearest stream. Once in the stream, existing flood 
routing algorithms are probably sufficient, at least for an 
initial approximation to the runoff distribution. In the 
coming field season we plan to use our arrays of thermistors 
and soil moisture blocks, plus thermal imagery from 
Landsat-C, HCMM (if available) , and NOAA environmental 
satellites to estimate water characteristics of the snowpack. 


r.ITBPRTUPS P’^VIEW oy REMOTE SENSING APPLICATIONS TO 
HYDROLOGIC MODELING 

1. Development of Bihliograpliy and Computerized 
t^etrioval system 

T!ius '^ar our efforts on this task have been devoted to 
gathering irticles on the specific problems that we have 
coht'tonted, and to developing a computerized information 
rotrievcvl system for this bibliography- In comparison to 
other literattire data bases, this one is quits small and 
rather specialized- At present it deals only with those 
aspects of hydrologic modeling where we have investigated the 
role of remctw sensing- 

The purpose of the literature retrieval system is 
primarily to organize the information in the library. The 
first task was to define a list of keyword topics into which 
the articles and books could be categorized. For each 
article, book, and proposal a card was made listing the 
author or authors, the title, the date of publication, the 
publication title for a journal or the publisher for a book, 
the topic or topics of the entry, and any other relevent 
information- The cards are filed alphabetically ^by author 
and kept as a complete catalogue of the library- They are 
computerized in order to facilitate their rapid sorting by 
author, topic/ or publication- 

CAHDCAT, the card catalogue program, is capable of 
maintaining a library of up to about 30,000 bibliographic 
entries. (At present only about 500 entries are supported.) 
It can add entries to this library and link them together in 
four different ways. Thus the user can obtain a listing of 
the entire data base alphabetically by author, of a given 
topic alphabetically by author, of a given author 
chronologically (latest first) , or of a selected entry by 
specifying a title. Entries may be deleted at any time. The 
user may also specify unions and intersections of topic 
listings so that subtopics may be used in the cataloguing 
process- For example, entries under subtopic “a" may be 
obtained. A listing of all entries under both topics "3" and 
”a" would contain those entries under topic "B," subtopic 
"a." By requesting a listing of all entries in the union of 
topics "A" and "B" the user would obtain the entire list of 
entries for both topics. Finally, the user may obtain a 
listing of all entries by a single author- When not 
executing, the CARDCAT program stores the data file on disk. 


2. Future Remote Sensing Applications to Hydrologic 
Modeling 

Current inputs to hydrologic modeling from remote 
sensing have mostly been limited to presence or absence of 
hydrologic phencmena- This is evident in the extensive use 
of remote sensing data for mapping snowcover areal extent 


NL. 
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(e.q. Barnes and Bowley, 1970; Siefert et al, 19‘75; Eango, 
1975), Improved accuracies of runoff prediction from the use 
of this application have been striking (Rango and Salomonson, 
197 5), 


However, most operational hydrologic models are lumped 
parameter models, usually based on one or more convienient 
indexes (such as snow course measurements). Remote sensing 
data offers high resolution, multispectral information that 
cannot be effectively used by such models. With the 
development of distributed paramenter hydrologic models, such 
as an energy balance snowmelt model, comes the opportunity to 
more fully utilize remotely sensed data. 

Techniques developed for remotely sensing thermal 
characteristics of the surface and the atmosphere have been 
restricted to oceanic regions due to noise problems from 
terrain and the differential thermal properties of land 
materials (Yates and Bandeen, 1975), It is hoped that these 
techniques can be extended over an alpine snowfield where the 
thermal properties of the surface are known and where the 
terrain effects can be filtered. As new remote sensing 
devices become operational, such as the HCMM and landsat-C, 
the effective use of thermal data inputs for hydrologic 
modeling will become more important (Lockheed Electronics 
Co. , 1977) , In this regard, we will be reviewing literature 
on the following topics; 

i. Remote determination of soil and snow 
moisture content, 

ii. Plant canopy temperature determination and 
measurement, 

ill. Determination of cloud cover thickness, 
temperature, and water content- 

iv. Techniques for thermal mapping of 
snowfields; especially mapping the isothermal (OOG) 
portions of a snowfield, 

V, Determination of atmospheric temperature 
and vapor pressure profiles. 
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III- SOFItfJ.SE DEVEIOPMEMT - 

NASA HYDROLOGY PF(CJ3CT PHOGHAH SISUCTURS 

The NASA hydrology project programs are designed to 
perform specific tasks, and to write the output in a file 
which may in turn be accessed by other programs. In order 
that the hassles involved be minimal, they must be 
standardized in output format. 

By agreement then, we adopt the following standards: 

1. With some exceptions, output files are in 
binary, and are written as sequential files on tape or 
disk. 

2. The first record should contain the dimensions 
of the output file (e.g,. NX by NY for a grid) and other 
appropriate information, such as latitude and longitude, 
grid spacing, etc, 

3. Gridded information is written by row, starting 
from the SW corner, Extensive existing software (e. g, 
DIRS) which does not use this convention should not be 
rewritten though. Instead, the output files can be 
subseguently re-formattad. 

All progiants should be briefly documented. The 
documentation should contain the input specifications, a 
description of the methods used, and an explicit description 
of the output records- Programs should also give a sample 
print-out of the output data set. Because we use a variety 
of languages, we must take care to make data sets compatible. 
Load modules for complete programs should be stored in the 
RDS HYDROLOGY. PROJECT, PROGLIB on D10300 unless you have a 
particular reason for putting them elsewhere- This disk 
should be otherwise used for data sets only, 

I, Topographic analysis, 

A- DTTGRG - selects data from digital terrain 
tapes and writes into a disk or tape file. Options for 
interpolaticn and smoothing- 

B. TBERAIN - calculates slope, exposure, horizons, 
latitudes, longitudes and elevations. Writes into file. 
Can be subgrid of that written by DTTORG, 

II, Solar radiation, 

A. RADPENT - prints out in table format 
measurements with Landsat radiometer. Writes results 
into seguential indexed file for input to other solar 
radiation programs. 
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B, SOLGRD {not yet complete) 
radiation. Numerous options: 


- calculates solar 


1. Selected points in selected wavelengths at 
selected time step- Print out- Used for 
verification and calibration. 


2- Selected wavelength intervals at specified 
time over entire grid- Used for comparison 
with imagery, 

3. Selected wavelength intervals at specified 
time step over entire grid. Used for model 
in p ut- 


il. Selected integrations over desired 
intervals, 


Options 3 and 4 may be in terms of net solar 
radiation as well as incoming. Program uses either 
estimated or measured reflectances; estimaxed or 
measured absorption and scattering coefficients and 
portion of radiation scattered forward. Interruption by 
canopy cover can be specified for all options. 

ill. Longwave radiation- 

A. LCNGP5NI - prints out in table format longwave 
radiation raeasureraentso 

B, SPACSAVR - (formats measurements for input to 
other programs) ; writes binary data to disk in packed 
format- 

C- TZIEMP - tirae/elevation interpolation to 
estimate Kumidity and temperature at all elevaMons and 
times, 

1, Calculates lapse rates for temperature and 
huraidity. 

2. Adds pressure/ temperature, and humidity 
measurements from Bishop weather station, 

3- plots time-space graphs for temperature 
and humidity. 

4, Writes interpolated temperature and 
humidity to disk for IZGRID, 

D, INRLUX - calculates longwave radiation from 
temperature, humidity, and topographic data. Uses 
measurements for calibrations. Options: 

1- Selected points at selected time step used 
for verification and calibration. 


2- At speciflad time aver entire grid. 

3. At specified time step over entire grid. 

Used for model input. ' 

4.. Selected integrations over desired time 
intervals over entire grid. Used for model 
input. 

Temperature and 'humidity data. 

A- TZGEID - coarse model. Uses TZT3MP to specify 
temperature and humidity values at every grid point for 
every time step. 

Hind information - not yet developed. 

Snowmelt modeling. 

A. SHOWBAl - (from A"nderson, 1976) uses input fro 
previous programs to run snowmelt model for all grid 
points. 

Satellite data analysis programs. 

A, DIRS - NASA program for rectification, 
rotation, and printer display of Landsat data- System 
elements: 

1, COMCCT - Converts rectified image or 
subimage from composite to CCT format- 

2. COMPOS - Refcrraats CCTs into composite 
tape format; produces shade prints; does edge 
image correlation- 

3. DTCOSYS (not yet implemented) - Accepts 
film recorder format image tapes and applies 
annotation contrast and expansion operations 
to the image; produces film recorder output 
tapes- 

4, HAPX (not yet implemented) - Extracts 
single band quadrilateral subimage areas from 
composite format rectified image tape, 

5. MAP2X - Rectifies quadrilateral subimage 
areas from coraposita tape with GSOM file and 
writes output in composite or film recorder 
format, 

6, QUIKIOOK (not yet implemented) - Accepts 
composite format image tapes and produces film 
recorder format output image tapes with 
optional geometric, radiometric, and 
annotation enhancements. 


-- 22 - 


7. RECTFI - Cross checjcs control points (GCPs 
and CCPs) ; computes global mapping functions; 
creates interpolation grid; writes GEOH file 
on composite tape; optionally rectifies full 
images 

8 - BSFRA^B - Vertically mosaics unrectified 
Landsat scenes which are contiguous. 

9- REGISTER (not yet inipleraented) - Converts 
correlation control points into pseudo GCPs by 
computing UTM coordinates at point and 
storing point in GCP format, 

10. SIATDUi’IP (not yet implemented) ~ Extracts 
image format center time and attitude estimate 
from SIAT file. 

B, RABTAL - reformats BIBS output and converts to 
radiance values, 

VIII, Mapping programs, for interpolation and display of 
results. 

A. BITMAP - maps output from DTTORG, as contour 
maps or perspective plots, 

B. NCAH graphics package (not yet implemented) - 
powerful set of subroutines for graphic output on 
Calc o rap plotter, 

C. NBWGRD - transforms irregularly spaced data 
into square or rectangular grid by using gravity model 
to form irregular grid and smoothing with Laplacian; 
bicubic spline interpolation then creates grid, 

IX, Miscellaneous subroutines, in addition to those in the 

IMSL package. These are imbedded in the above programs, 
but others may find them useful separately, 

AHGIN - solar geometry for slope 

ASPECT - slope and exposure for point on grid 

COF 4 FIG - calculates absorption and scattering 
coefficients and portion of radiation 
scattered forward from set of three global 
radiation measurements 

COQISS ~ Fast Fourier Transform 

CSTEP - Fourier Transform of a step function 

EPHEM - declination, radius vector, and 

equation of time as function of date 


ilORTAE ** calculates horizon angle file for 
point in grid 

IDMSS - converts degrees, minutes, seconds (as 
single integer) to seconds 

MOViE - Moves one vector into another (in 
Assembler - very fast) 

NEHTOM - root finding by Newton's method 

RSINTG ^ Rayleigh scattering coefficient 
averaged over specified wavelength range 

F.SIN'TE “ Rayleigh scattering coefficient at 
specified wavelengths 

SATN and SA"I - saturation vapor pressure over 
water and ice 

SCXN1G - integrated solar constant in 
specified wavelength range 

SCI NT F ~ solar constant at specified 
wavelengt hs 

SKYFAC - diffuse solar radiation and 
atmospheric thermal radiation view factor 

SNPATH “ sun position for date, latitude, and 
time 


ZFEO - zeroes a vector (in Assembler ” very 
fast) 

ZPSESS - calculates pressure as a function of 
elevation, temperature, and lapse rate 


XV. ACGOHPLISnaENI’S FOHE!S:]SH B'i END OF FIRST Y3AR CF GRANT 
AND RFCOMHENDATIOMS FOR FOETHER WORK 

Prior to the end of the first year of the grant, we 
foresee completion of the following tasks: 

1. The clear-sky solar radiation model should be 
developed, with the inclusion of reflectance from adjacent 
slopes. The model will be adeguately tested. 

2. A method of determining snow albedo from Landsat 
data will be completed and tested. 

3- A method of calculating clear-sky longwave radiation 
at varying altitudes will be completed and tested. 

4- A geocoded hydrologic information system will be 
developed for the study area. 


Further tasks which we hope to accomplish in the next 
two field seasons are; 

1, Interpolation of Soviet satellite imagery of the 
Fergana Valley test site, 

2, Review of hydrologic modeling literature, to 
evaluate what model changes are needed for remote sensing to 
play a more meaningful role. 

3, Development and testing of frost penetration 
calculations for the study area- 

4,, Development of algorithms to couple snowmelt model 
to flow model, 

5. Use of thermal imagery from Landsat-C,HCHH, and 
NOAA environmental satellites such as TIROS N to better 
estimate longwave radiation from clouds and forests. 

6, Use of thermal imagery to estimate wind distribution 
over rugged terrain. 
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